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ABSTRACT 

The First Stars in the Universe form out of pristine primordial gas clouds 
that have been radiatively cooled to a few hundreds of degrees Kelvin either via 
molecular or atomic (Lyman-a) hydrogen lines. This primordial mode of star 
formation is eventually quenched once radiative and/or chemical (metal enrich- 
ment) feedbacks mark the transition to Population II stars. In this paper we 
present a model for the formation rate of Population III stars based on Press- 
Schechter modeling coupled with analytical recipes for gas cooling and radiative 
feedback. Our model also includes a novel treatment for metal pollution based 
on self-enrichment due to a previous episode of Population III star formation in 
progenitor halos. With this model we derive the star formation history of Popula- 
tion III stars, their contribution to the re-ionization of the Universe and the time 
of the transition from Population III star formation in minihalos (M ~ 10 6 M Q , 
cooled via molecular hydrogen) to that in more massive halos (M >2x 1O 7 M , 
where atomic hydrogen cooling is also possible). We consider a grid of models 
highlighting the impact of varying the values for the free parameters used, such as 
star formation and feedback efficiency The most critical factor is the assumption 
that only one Population III star is formed in a halo. In this scenario, metal free 
stars contribute only to a minor fraction of the total number of photons required 
to re-ionize the universe. In addition, metal free star formation is primarily lo- 
cated in minihalos and chemically enriched halos become the dominant locus of 
star formation very early in the life of the Universe — at redshift z ~ 25 — even 
assuming a modest fraction (0.5%) of enriched gas converted in stars. If instead 
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multiple metal free stars are allowed to form out of a single halo, then there is 
an overall boost of Population III star formation, with a consequent significant 
contribution to the re- ionizing radiation budget. In addition, the bulk of metal 
free stars are produced in halos with M > 2 x 10 7 M Q . 

Subject headings: cosmology: theory - galaxies: high-redshift - early universe - 
ISM: evolution - stars: formation 



Introduction 



Population III stars are considered to be the first luminous objects formed during the 
Dark Ages of the Universe, when the hydrogen is in a neutral state (e.g. see lBromm fc Larson 
2004 ). The first generation of stars, formed out of pristine primordial gas, had a top-heavy 
initial mass function, with a typi cal mass scale o f order of lOOM^ and most probably 
just one star per halo (e.g. see lAbel et al.l 120021 ; lO'Shea fc Normanl 120071 ). These stars 
start forming afte r about 30 — 40 million years fro m the big-bang a t redshift z ~ 55 — 65 
(INaoz et al.l 120061 ; iTrenti fc Stiavellil 120071 : see also iGao et al.l 120051 ) and, given their high 
mass, they live only a few million ye ars ending with ei ther a pair instability supernova phase 
or a direct collapse to a black hole (IHeger et al.l 120031 ). 



Population III stars thus initiate the chemical enrichment of the Univ erse and open the 



way to more normal modes of st ar formation, namely Population II (e.g. see lOstriker &: Gnedin 
19961 : iFurlanetto fc Loebl 120031 ) . In fact, the metals released into the IGM after a pair in- 
stability supernova explosion can trav el outside the parent dark matter halo that hosts the 
Population III star. Calculations by iBromm et al.l (120011 ) found that a region containing 
up to about 1O 8 M can be enriched to a critical metallicity Z crit > 10~ 4 Z Q by the most 
massive pair instability supernovae. More typical explosions may instea d enrich significantly 
less gas (~ 10 6 M^) although at a correspond i ngly higher metallic ity (see lBromm et al.ll2003l : 



Kitayama fc Yoshidall2005l : iGreif et al.ll2007l : IWhalen et al.ll2008l ). Even in the latter case, a 
halo of mass < 1O 8 M that had one of its progenitors hosting a Pair Instability Supernova is 
still like ly to be enriched t o an average metallicity of > 10~ 4 Z Q thanks to violent relaxation 
mixing (ILynden-Belll Il967l ) during its hierarchical build-up. 



Population III stars are also the sources that start to re-ionize the Universe, creating 
ionized islands within the neutral hydrogen inter-stellar and inter-galactic medium. Ionizing 
photons are emitted with an enhanced efficiency compar ed to Population II stars due to the 
high effective temperatures of massive metal-free stars (ITumlinson fc Shulll l2000l : ISchaerer 
20021 ) . and these sources could be responsible for a significant fraction of the Thompson 
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optical depth to reionization deriving from z > 7 (jShull & Venkatesanll2008l ). Another hint 
suggesting that Population III stars contribute significantly to the re-ionization of hydrogen 
can also be inferred by the rapid evolution of the galaxy luminosity function at z > 6, which 
im plies that observed galaxies alone do not seem capable of re-ionizing the Universe (e.g. 



sec 



Oesch et al.ll2008l : Bouwens fc Illingworthll2006r ). 



Two main modes of Population III star formation have been proposed: either in mini- 
halos with virial temperature of T vir rs 10 3 K, where the gas is cooled via molecular hydrogen 
{H 2 ), or in more massive, rarer, halos with T vir lO 4 ^, where cooling t hrough atomic hydro- 



gen (Lyman-a) lines becomes possible (e.g. see lBromm &: Larsonll2004l ). H2 is formed during 



the initial collapse of the gas within the minihalo, but it is sensitive to photo-dissociating 
radiation in the Lyman Werner band ([11.18 : 13.60]eV). Thus, in presence of a sufficiently 



Ciardi et al. 2000; Glover fc Brand 2001; Machacek et al. 2001 
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Interestingly, the main coolant of halos with a virial temperature of T V i r « 10 4 K in a strong 
LW bac kground may continue to remain molecular hydrogen according to numerical simu- 
lations: I O'Shea fc Norman! (120081 ) found in fact that while in presence of a strong feedback 
T vir « 10 A K is required for collapse, still the cooling is driven at the c enter of the halo b y 
molecular hydrogen, formed thanks to the high central density (see also lWise & Abelll2007l ). 
The feedback induced by nearby Population III sourc es may also be positive, as, e.g., a 



soft X -ray background enhances the H 2 production rate (IRicotti et al.ll2001t iMachacek et al. 



20031 ). Therefore the termination of the first epoch of Population III star formation depends 
critically on the relative weight of these two competing process, which in turns is influenced 
by the local topology of the IGM, by the spatial distribution of the sources and by their 
IMF (which affects the relative efficiency of X ray to Lyman Werner photon production). 
In addition to radiative feedback, Population III star formation is influenced by chemical 
feedback. This can be broadly classified as (i) self-enrichment due to a previous episode of 
star formation in a progenitor of the halo considered and (ii) metal pollution due to galactic 
winds originated in a nearby halo. Given such a complex scenario it is not surprising that in 
the literature there have been many investigations focused on characterizing the Population 
III star formation rate and the nature of the transition from Population III star formed in 
minihalos t o Population III stars formed in T,„> ~ 10 4 halos and from Population HI to Pop 



ulation II flMackev et all 12003 : iFurlanetto & Loebll2005t iGreif & Bromml 12009 : ISmith et al. 
2008h . 



The formation of Population III stars is typically investigated by means of two comple- 
mentary approaches: (i) analytic models aimed at deriving an average s tar formation rate 
- these usually rely on a dark matter halo formation rate derived with a lPress fc Schechter 
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( 119741 ) lik e formalism combined with recipes to populate the dark halos with Population 
III stars jMackev et aljl2003l : iGreif fe Bromml hood : IWvithe fc Chenll2006h : (ii) high reso- 
lution hydrodynamic - radiative transfer simulations that follow in detail the collapse and 



the early stages of formation of a single Populati on III star (e.g. see lAbel et al.l 12002 



O'Shea & Norman! I2007L 120081 : lYoshida et al. 1120081 ). 



In this paper we have two main objectives. First, we focus on the characterization of the 
global star formation rate of Population III stars, thus adopting an analytical approach. We 
resort to physically motivated recipes to identify the conditions under which it is expected 
that the primordial gas within a dark matter halo can cool and trigger a gravitational 
instability which leads to a protostellar core. These recipes include the effects of a photo- 
dissociating Lyman Werner background derived both self-consistently from our model as well 
as by adopting a reference reionization history of the Universe. 

Our second goal is to quantify the probability that a newly formed dark matter halo 
with virial temperature T V i r ~ lO 4 ^ has been previously enriched by one or more episodes 
of H 2 Population III formation in one of its parent halos. Chemical enrichment of such halos 
is in fact crucial not only to assert the relative contribution of Population III star formation 
via the ato mic and molecular co oling channels, but also to evaluate the formation rate of 
quasistars (IBegelman et al.l 120061 ) at z ~ 15 which have been proposed as progenitors of 
the supermassive black holes present after the end of rei onization. Quasista rs are in fact 
able to form only if the gas is not polluted by metals (lOmukai et al.l 120081 ) . Our novel 
approach to self-enrichment is based on the properties of the Gaussian random field of the 
primordial density fluctuations, which allow us to derive a closed form for the probability 
that a dark matter halo of mass Mi at redshift z\ had at redshift z 2 > z\ a progenitor of mass 
M 2 > M\ (ITrenti fc Stiavellil 120071 ). We the n combine our resul t s on self-enrichment with 
the probability of wind pollution derived by iFurlanetto fc Loebl ( 120031 ) to infer the overall 
likeliness of collapse of pristine gas in halos with T vir pa 10 4 i^. 

This paper is organized as follows. In Sec. [2] we introduce our model for Population 
III star formation, including radiative and self-enrichment feedback; The model is applied 
in Sec. El to derive the global Population III star formation rate and in Sec. H]to obtain the 
enrichment probability of T vir pa 10 A K dark matter halos. Sec. O discusses the implications 
in terms of contributions to reionization from Population III stars and Sec. [6] concludes. 
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2. Population III star formation model 



To derive the star formation rate of Population III stars we combine the dark matter 
halo formation rate with an analytical model to populate halos with stars. In this paper 
we assume a flat concordance AC DM cosm ology, with the cosm ological parameters given 
by the WMAP Yr5 best fitting parameters iKomatsu et al.ll2008h : tt A = 0.72, Vt m = 0.28, 
Qb = 0.0462, as = .817, n, s = 0.96, h = 0.7. We also assume a primordial helium mass 
fraction Y = 0.2477 JPeimbert et~al1l2006h . 



2.1. Minimum Minihalo Mass for Population III formation 



The minimum dark matter halo mass M capable of cooling by molecular hydrogen at 
redshift z is estimated by re quiring the cooling time r coo ;(M, z) to be no larger than the local 
Hubble time tn(z) (e.g. see lCouchman fc Reeall986h . We write the cooling time as: 

/w \ 3k B T vir (M,z) 

T cool {M,z) = — — — — , (1) 

2A(T vir ,n H )f H2 

where ks is Boltzmann's constant and T V i r (M,z) is the virial temperature of the halo, A 
is the cooling function per Hi molecule, which depends on the temperature and on the 
hydrogen number density nn and fn^ is the molecular to atomic hydrogen fraction. We 
write T vir following iTegmark et al.l (119971 ) as: 

2/3 



T vir (M,z) ~ 255AK 



M 



(l0 6 M G 



l + z 
31 



(2) 



For th e molecular hydrogen cooling function A, we use the form derived by I G alii &: Palla 
( 119981 ). which we approximate between the temperatures of 120 K and 6400K (the range we 
are interested in for cooling in minihalos) with: 



A(T,n H ) ~ 10 



-31.6 ( JI_Y A x 
\100K J 



10" 4 cm- 3 



ergs 



(3) 



We estimate the hyd rogen number density in a halo from its virial density (e.g. Eq. 22 of 
Tegmark et al.lll997l ) to find: 

3 



n H ~ 1.01 



l + z 
31 
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Replacing Eqs. El El H in Eq. [TJ we find: 



t coo1 ~ 3.46 x 10 10 s 



/ M 
V10 6 M C 



cm 



-1.6 /, \ -5.4 

l + z 
31 



-i 

H 2 - 



(4) 



(5) 
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We can then obtain the molecular hydrogen fraction required for cooling by equating the 
cooling time given by Eq. [5] to the local Hubble time approximated at z ^> 1 as: 

, , 1.191 x 10 15 s (\ + z\~ 3/2 nn „ - fl + zY 3/2 

hJUZ hr - 3.216 xlO-s — . (6) 



/ M \ Lfl (\ ■ : x ' 



This gives us: 

Tegmark et al.l (119971 ) determine the (maximum) molecular hydrogen fraction capable of 
forming in a halo as: 

/ rp \ 1.52 

( IbW J • (8) 

Equating the required molecular hydrogen fraction for cooling within a Hubble time given 
by Eq. [7] with the maximum that can form (assuming T = T vir ) we find the minimum mass 
for a minihalo in order to cool within a Hubble time {M tH __ coo i), namely: 

/i , \ -2.074 

M tH . cool ~ 1.54 x 1O 5 M ( J . (9) 



/fl2,mttx ^ 3.5 X 10" 



2.2. Minimum iJ 2 cooling mass in presence of radiative feedback 

In order to compute the effect of a radiative flux in the Lyman- Werner band on the 
formation rate of molecular hydrogen an d on the cooling of a minihalo, we resort to an 



approach based on iMachacek et al.l (120031 ) . We obtain the minimum halo mass capable of 
cooling via H 2 in the presence of a LW background by equating the timescale for photo- 
dissociation of molecular hydrogen (r diss ) to its formation timescale (rj orm ). 

In presence of a LW flux Flw = 47rJ 2 ilO~ 21 ergs~ 1 cm -2 Hzr 1 (whose c alculation is pre 



sented below in Sec. 12 . T[) . the dissociation timescale can be written as (IMachacek et al. 



200J): 

7.16 x 10 n s . . 

Tdiss ^ j , (10) 

The H2 formation time scale is given by: 

form k 7nHne k-jTL e ' 

where k 7 w 1.8 x io~ 18 T 88 cm 3 s _1 is the H~ formation rate which dominates the formation 
of molecular hydrogen and n e is the electron density. n e is obtained assuming a residual 
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ionizing fraction 2 x 10 4 ( [Peebles! 119931 ). Imposing the equilibrium molecular hydrogen 
function to be the minimum needed for collapse as given by Eq. [7] we finally find: 



M: 



H2—C00I 



6.44 x 10 6 MqJ^ 457 



l + z 
31 



-3.557 



(12) 



We note that Eq. [12] is in good agreement with the results by iMachacek et al.l (120031 ) at 
z ~ 30. Our derivation does however include an explicit redshift dependence . Eq. H21 also 



compares well with the results of numerical simulations by I O' Shea fc Norman! (120081 ). which 
include LW background of varying intensities. The redshift dependence w hich we find in- 
crease s the minimum mass required for cooling with respect to the formula by lMachacek et al. 
( 120031 ) at z < 30 and th is gives us a better agree ment with the numerical results, obtained 
for z < 25 (see Fig. 3 in lo 'Shea fc NormarjEjosh . 



In conclusion, for a dark matter halo to be able to cool via H2, its mass must be above 
both the limits set by Eq. [9] and Eq. [121 that is 



M„ 



max(M tH -cooi; M 2 



Hi— cool ) 



(13) 



2.3. 



Cooling in halos with T vir > lO^K 



Pristine halos with a yirial temperature above T vir ~ 10 4 i^ can cool irrespective of the 
LW background intensity (jO'Shea fc Normanll2008l ). In fact, not only ato mic hydrogen cool- 
i ng be comes in principle available, but also cosmological simulations by lO'Shea fc Norman 
(120081 ) have shown that a small fraction of H2 can still be produced at the center of the halo 
thanks to the high density and self-shielding of the surrounding gas. Once the gas temper- 
ature starts to decrease, further cooling and collapse will proceed progressively faster via 
molecular hydrogen as the halo temperature is initially high enou gh to enhance the abun- 
dance of H~, a precursor for H 2 production (ILepp fc Shulll Il984h . In our model we thus 
consider that all halos above the Mt=w 4 k limit will cool efficiently: 



M T=10 4 K = 7.75 • 10 6 M 







l + z 
31 



-3/2 



(14) 



2.4. Forming Population III stars 

There is of course a delay between the virialization of a dark matter halo potentially 
able to cool via H2 — that is of mass M > M m i n (Eq. [T3l) — and the actual formation of 
a Population III star. We estimate this delay considering two contributions: (i) the actual 
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time needed to cool down to a few hundreds degrees Kelvin and (ii) the free fall time for the 
gravitational collapse once cooling has triggered the Jeans instability. 



The H 2 cooling time can be obtained by Eqs. and [SJ 

-2.627 

-,„„/ ■ 2.3S • lO'-'sf 



1O 6 M / V 31 



-6.94 



(15) 



while the free fall time can be obtained from the Jeans instability timescale, taking into 
account that during the cooling phase the density of the gas has increased by about a factor 
7: 



tf f = 2.77 x 10 14 s 



l + z 
31 



-3/2 



(16) 



Therefore a Population III star originated in dark halo virialized at redshift z v i r will be 
formed at redshift Zf orm < z V i r such that: 



tcooi + tff = 3.21 x 10 15 s 



1 + Z fo 

31 



-3/2 



1 + Z V 

31 



-3/2' 



;i7) 



where the right side of the equation simply derives from the age of the Universe at z ^> 1 
(Eq.[6Tj. 



2.5. Metal enrichment probability 

In order to account for previous episodes of Population III star formati on in a progenitor 



of a h alo of mass M at redshift z±, we resort to the method presented in iTrenti fc Stiavelli 



(120071 ). based on the linear growth of density perturbation in the context of spherical collapse. 
We start by assuming that a newly virialized halo has an average linear overdensity 5(z) = 
1.69 as estimated by a top-hat filter on a scale R = (3M/4vr(p)) 1 / 3 . Then for a progenitor 
mass M prog < M we compute the extra variance in the density power spectrum a 2 add = 
a 2 (M prog ) — a 2 (M) and then the refinement factor N re f = M/M prog . With these ingredients 
we can write the probability distribution for the maximum of N re f Gaussian random numbers 
with variance a^ dd as the derivative of the N re f power of the Partition function for a normal 
distribution with zero mean and variance o\ dd . In the context of spherical collapse this 
translates to a probability distribution for the formation redshift of the first progenitor of 
mass M prog of a halo M virialized at z\. 

For every progenitor mass M prog we then compute the delay time (tff + t coo i) needed to 
form a Population III star in the parent halo and from this we derive the minimum redshift 
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{zmin.seed) at which such a halo must form in order to pre-seed the descendant halo. Of 
course, for some values of M prog the delay time might be longer than the Hubble time, this 
simply means that no pre-seeding is possible from parent halos of mass M < M prog . We then 
integrate the probability distribution for the formation time of the parent of mass M prog for 
z > Zminseed to obtain the preseeding probability from a progenitor at this mass scale. The 
overall probability of preseeding is the maximum preseeding probability computed over all 
the possible progenitor masses. 



2.6. Prom dark matter to stars 



The dark matter halo formation rate is derived in our reference model using the lSheth fc Tormen 

( 1999 ) mass function. The Sheth fc Tormen ( 1999h m ass function is in better a greement with 

Press fc Schechterj ( 1974 ) mass function at z < 30 ( Heitmann et al. 



N-body simulations th an the 



20061 : iReed et al.ll2006l ). Note that differences of ~ 20% have been observe d between the mea- 



suremen ts from cosmological simulations and the ISheth fc Tormen! (119991 ) mass function and 



that the 



Wa rren et al.l (120061 ) ma ss function appears a be tter matc h to the numer i cal re sults 



( ILukic et all 120071 ). However the ISheth fc Tormenl (ll999l ) and the IWarren et all (120061 ) for- 
mulae give very similar results in the r ange of halo masse s of interest for P opulation HI star 
forma tion (M < 1O 8 M - see Fig. 3 in lLukic et al.ll2007l ) thus we keep the lSheth &: Tormenl 
( 119991 ) model as our reference. We then compute the formation rate of H% Population III 
stars by integrating between M min (z) and M T=W 4 K (z) the number of halos per unit mass per 
unit redshift dN(M, z) / dMdz, convolved with the probability that such halos are pristine 
(see Sec. 12. 5p . 

The characteristic mass of Population III stars and the form of their initial mass function 
are highly uncertain, even though they are likely very massive — of the order of O(1OOM ) 



(e.g. see 



Bromm &: Larsonll2004t). This expectation is based on theoretical models and nu- 



meric al simulations (lAbel et al.ll2002l ; lOmukai fc Yoshiill2003l : lYoshida et al.ll2006l ; iGao et al. 



20071 ). but there is some tension with the abundance patterns observed in the most metal 
poor Milky Way stars, which are better explained under the assumption that th eir progeni- 
tor P opulation III stars were only moderately massive — 8M Q < M < 42M — ( iTumlinson 



20061 ). There i s however no guar antee that the progenitors of the extremely metal poor stars 
considered by iTumlinsonl (120061 ) are Population III stars formed before the reionization of 
the Universe: if their progenitors formed instead in presence of a strong UV background 
within a reionized bubble, then thei r expected mass is a bout ~ 40M Q fully consistent with 
the inference from the observations (lYoshida et al.l 120071 ). Within this complex scenario we 
choose to adopt conventionally one Population III per minihalo (jO'Shea fc Normanl 120071 ) 
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and we consider a ISalpeterl (119551 ) mass function in the range [50 : 300] M (average mass 
1OOM ), as suggested by the theoretical investigations. A modification in the initial mass 
function used in our model primarily affects the enrichment history of the IGM and thus the 
transition toward Population II star formation. If Population III stars are less massive than 
we assume, then the efficiency of metal pollution may b e reduced as core c ollapse supernovae 
explosions are not as energetic as pair instability ones (IHeger et al.l 120031 ). A mass function 
more biased toward very massive stars with M > 27OM would also reduce the efficiency of 
metal pollution, because these stars directly collapse into black holes without an explosive 
phase (IHeger et al.l 120031 ) . 



The formation rate for Population III stars in halos with T > 10 4 K is similarly computed 
from the dark matter halo mass function for M > Mt =W 4k, again after convolution with 
the probability that the gas forming new halos has not been contaminated by metals. As 
no sign of fragme ntation has been found du ring the collapse of metal free halos with masses 
up to 2 x 10 7 M Q (jO'Shea fc Normanl 120081 ). we adopt the same initial mass function as for 
H 2 cooled Population Ills (one star per halo, Salpeter in [50 : 3OO]M ). The efficiency 
of star formation is however uncertain and thus we explore different models where multiple 
Population III stars in a single halo are allowed, adopting a reference star formation efficiency 
(star to gas mass ratio) of ep op ni = 0.005 and epopin = 0.05. 

Finally the star formation rate in enriched gas (Population II stars) is computed by 
convolving the dark matter halo mass function with the preseeding probability (the comple- 
mentary of the pristine probability) and assuming a star formation efficiency of ep op u = 0.005 
or epopii = 0.05 (to explore the uncertainties in this parameter), a Salpeter mass function 
from in the range [1 : 100]M o (average mass ~ 3M ) and an average metallicity 10~ 4 Z Q . 
Our choice of the IMF for metal enrich ed gas reflects th e expectation that the typical mass 
of stars was higher at higher redshift (iTumlinsonl 120071 ) . In our model, the IMF of metal 
enriched star formation impacts the radiative LW feedback, but only in a minor way due to 
its self- regulating nature (see Sec. [3]). 



2.7. Flux in the LW band 

We compute the LW flux that enters in Eqs. [10] by means of two different approaches: 
(i) self-consistently from our model, based on the star formation rate and (ii) adding a 
reference number of LW photons to those derived self-consistently in order to take into 
account Population II formation not included in our model. The first method is most suitable 
at z > 20, when Population III are most likely the dominant sources of radiation. At 
lower redshift protogalaxies become more and more common, the Universe starts becoming 
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reionized and our simple model for star formation does not capture all the Population II 
formation that is available, therefore using a reference LW photons production provides a 
good check on the validity of our assumed LW flux. 

For the self-consistent LW flux calculation, we obtain that Popul ation HI stars f ollowing; 



our assumed IMF emit a LW flux that is 7.5% of the ionizing flux (jSchaereri 120031 ) . Metal 
enriched stars have instead a higher ration of LW to ionizing photons, because these stars 
have a lower effective temperature, but their LW photon yield per uni t mass is also lower. 



We assume the following photon yields over the star lifetime (based on ISchaererl 120031 ) : 

1. Population III: 8 x 10 60 Mq 1 ; 

2. Population II: 8 x lO^M" 1 ; 

The comoving LW photon density tilw is then associated to a flux: 

J 21 = 1.6 x 10- 65 ( -^—] ( —) 3 ergs^cm^Hz-V- 1 . (18) 
yiMpc-i J V 31 J & v ; 

In order to compute tilw( z ) we only consider star- formation that has happened within 
a redshift interval such that the photons have not been redshifted out of the LW band 
on average. This means that the upper limit, expressed in term of the redshift is z up = 
12.39/11.18(1 + z) — 1. We also take into account the screen provided by primordial Hi 
present in the IGM ou tside virialized halos, which can absorbs LW photons. Following 



Trenti fc Stiavellil (120071 ) a flux of J 2 i = 1-58 x 10 -3 is needed to photo-dissociate a molecular 
hydrogen density of 10 -6 times the neutral hydrogen density. Of course, a (very) small 
residual fraction of Hi will still be present but the re-formation of H 2 outside virialized 
structures is strongly suppressed not only by the radiative feedback, but also by the decreased 
average density of the univer se compared to that at the time of primordial H 2 formation 



(z « 200 — see |Peebleslll993h 



For the LW flux based on a fixed reionization history of the Universe, we assume that 
the number of LW photons is that produced self-consistently by Population III plus a con- 
tribution from other sources which reaches n^w = 7 x 10 66 at z = 10 (corresponding to 
J21 ~ 4.9). As a model for the redshift dependence of n^w we take inspiration from the 
rapid growth of the fraction of mass in collapsed dark matter halos that can host stars and 
we write: 

n LW (z) = 7 x 10 66 x i 3-3-3.3x(i +2) /n (1Q) 

The flux in the LW band approximately matches the reionizing flux given for a st ellar pop- 
ulation with a Salpeter IMF in the range [1 : 1OO]M and metallicity Z « 1O~ 3 Z (ISchaererl 
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20031 ). Thus Eq. [19] implies that about one ionizing photon per hydrogen atom is emitted 



at redshift 10, a budget still well short of what is required for reionization after considering 
the dumpiness of the IGM and recombination. 



3. Population III star formation rates 

From our fiducial model, which has one first star per halo (see Tabled] for a summary of 
the main parameters), it is immediate to note that Population III stars initiate the chemical 
enrichment of the Universe well before redshift z = 50 (see Fig. []]). This result derives from 
the relatively small mass required at very high redshift to be able to cool via molecular 
hydrogen. Such mass is in fact as low as M ~ 4 x 10 4 M Q at z = 60. At lower redshift the 
minimum mass for cooling progressively increases but halos capable of cooling become more 
and more common, so the Population III star formation rate steadily increases. Eventually 
- around redshift z ~ 35 — radiative feedback in the LW band starts to significantly 
increase the mass required for cooling and the star formation rate of this class of Population 
III stars levels off at w lO~ 5 M Mpc _3 yr _1 , that is about one star in a comoving Mpc 3 formed 
per unit redshift. In our reference model the self-shielding mass becomes larger than that 
of a T vir = lO^K halo at z ~ 13. Following Sec. 12.31 we assume that such halos can cool 
independently of the LW background intensity. 

Forming Population III stars in these more massive halos is possible only if there were 
no previous episode of star formation within their progenitors. Thus, as long as cooling can 
be efficiently achieved via H 2 , halos with T vir > 10 4 K are most likely chemically enriched (see 
bottom right panel of Fig. [T]) and Population III stars in these halos are very rare compared 
to their counterparts in minihalos (see the upper left panel of the same figure). Their star 
formation rate becomes higher than that in minihalos only at z < 14, when the LW feedback 
strongly suppresses Hi cooling and thus it is more likely that a halo made entirely of pristine 
gas is able to grow via mergers to reach T vir > 10 4 before having the possibility to cool. 

The majority of halos with T vir > 10 4 K at z > 16 have been instead chemically enriched 
and Population II star formation grows rapidly in time. By z ~ 26 it becomes the dominant 
factor in the global star formation rate, despite our conservative assumption that only 0.5% 
of the gas is converted into stars (but note that this is still larger than the efficiency in 
Population III stars at z < 30 if only one per halo is formed). This early rise of Population 
II is a novel result which derives from our detailed treatment of self-enrichment. When 
metal transport is instead modeled via wind s, then transition toward Popu lation II stars is 



predicted significantly later, at 10 < z < 20 flFurlanetto fc Loebl 12003 . l2005f ) 
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The qualitative picture described in our standard model holds even if we vary the free 
parameters that regulates radiative feedback. In Fig. [2] we show the star formation history 
when LW feedback is strongly suppressed. Population III stars are formed at a higher rate 
at z < 30 compared to the standard model, but the difference is only a factor of a few (here 
the formation rate goes up to ~ 5 x 10 _5 M Q Mpc _3 yr _1 . This is about one fourth of the 
peak formation rate without any radiative feedback (see Fig. [3]). The main difference from 
the standard scenario is that Population III stars formed in halos with T vir > 10 4 K are now 
suppressed even at z < 15, because the enrichment probability of massive halos remains high 
(see bottom right panel in fig. [2]). The formation of Population II stars is instead essentially 
as in the standard model. Note that the picture does not deviate significantly from the 
standard model if we vary the efficiency of star formation in Population II stars, except 
of course for a corresponding proportional variation of their star formation rate. In fact 
Population III stars in minihalos are the main agents of the radiative feedback that leads to 
their suppression at z > 20 (see bottom left panel in Figs. [UTS]). Therefore their formation 
is self-regulated and tends to reach an equilibrium level. Similar results hold even when we 
add a stronger radiative feedback based on the fixed Lyman Werner background given by 
Eq. [H] (see Fig. @J. Note that in this case the background radiation greatly exceeds that 
created by Population III stars at z < 20 and thus the suppression of H2 cooling is even 
sharper. 



Fig. [5] shows a model with our sta ndard parameters except for the use of the lPress fc Schechter 



( 119741 ) mass function rather than the ISheth fc Tormenl (119991 ) formula. The predictions for 



the Population III formation rate are very different at z > 40, but once the self- regulated 
feedback phase starts the two models converge together. The strong difference at z > 40 
originates fro m the fact that Pop u lation III stars are hosted in very rare peaks at such early 



times: in the iPress fc Schechterl ( 119741 ) formalism these halos have v = 5^/a 2 (M) > 30 
and the ratio of the Sheth-Tormen to Press-Schechter mass function is proportional to 
exp(— 0.707z/)/ exp(— v) in the limit of very large v. An interesting open question is the 
form of the mass function for s uch rare peaks, which are expected to be progressively more 



spherical as they become rarer (IBardeen et al.l 119861 ). Fortunately the difference in the star 



formation rate does not propagate significantly below z < 30. 

A major qualitative change in the star formation history during the Dark Ages arises if 
we allow for multiple Population III stars in a single halo (see Fig. [6]). In this case the metal 
free star formation rate is comparable to that of second generation of stars (Population II) 
down to 2 ~ 10. From the star formation rate we can identify two eras: an early period 
dominated by Population III in minihalos up to z ~ 18 and a later period where larger halos 
are still able to form metal free stars. This happens because by allowing multiple Population 
III stars in minihalos their formation rate is significantly enhanced at later redshift over the 
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assumption of a single star per halo. In fact, when a significant LW background is present, a 
single minihalo can have up to 1O 6 M of gas, thus a ep ov m = 5 x 1CT 3 corresponds to a SFR 
50 times higher than the one obtained for a single metal free star per halo. In this scenario 
metal free stars formed in T vir = lO^K halos appear to be several orders of magnitude 
more common than in our standard scenario, again because the star formation efficiency is 
increased by more than two orders of magnitude compared to our standard model. They are 
therefore expected to dominate the production of ionizing photons (see Sec. [5]) until their 
formation is eventually expected to be terminated by chemical enrichment due to winds at 



4. Consequences for the chemical enrichment of T vir = 10 4 K halos 



The first generation of Population III stars — formed in minihalos — releases metals 
into the IGM and opens the way to metal enriched star formation when the gas recollapses 
later as part of a larger halo. When the second generation halo has a virial temperature 
T V i r w 10 4 K, the metallicity of its gas is expected to have been enriched to Z > 1O~ 4 Z , high 
enough to mark a transition toward a different mode o f star formation, especially if some 
dust is present ( iSchneider et al.ll2006l : Park et al.ll2008l ). Our model allows us to quantify 
the likelihood of this self-enrichment scenario as a function of the redshift of formation of a 
T V i r = 10 4 K halo (see bottom right panel in Fig. [1]). Interestingly, the probability of having 



a pristine halo large enough to reach T V i r = 10 4 K is small at very high redshift 



30) 



and progressively increases as the redshift decreases. This apparently surprising result can 
be understood in terms of the difference in halo mass required to cool via molecular and 
atomic hydrogen. At z > 30, the difference in the two masses is large (see upper right panel 
in Fig. [T]), hence it is highly likely that in the merger tree of the T vir = 10 4 K halo there has 
been a progenitor halo that was able to form a Population III star via H 2 cooling at an earlier 
redshift. However, as the redshift decreases, the LW self-shielding mass grows and it becomes 
progressively more unlikely that one of the progenitor halos hosted a star. Therefore halos 
with T vir « 10 4 K are more likely to be metal free at lower redshift. Our model does however 
not include metal pollution by winds, which becomes p rogressively more importa nt as the 



ripe 

redshift decreases. In fact, the wind pollution model by iFurlanetto fc Loebl (120051 ) predicts 



a sharp drop in the probability of forming pr istine halos around z ~ 15. Therefore if we 



combine our reference scenario results with the IFurlanetto fc Loebl (120051 ) model for the wind 
enrichment (e.g. see their Fig. 2), the overall picture is that the pollution probability for 
such halos remains high at all redshifts, probably presenting a minimum (with enrichment 
probability down to ~ 50%) between z = 15 and z = 20. This further strengthens our 
conclusion that Population III star formation in halos with T vir > 10 4 f^ is expected to be 
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subdominant compared to formation in minihalos if primordial stars are formed in isolation. 
This conclusion becomes stronger as the LW feedback efficiency is decreased or set to zero 
(see figs. [2] and [3]). 

A different possibility for star formation in massive metal free halos is the creation of a 
quasi-star, that is a bla ck hole formed via direct collapse accreting inside a massive envelope 



(IBegelman et al.ll2008l ). This scenario has been proposed to explain the rapid growth of 



supermass ive black holes, but a critical requirement for its viability is the absence of metals 



in the gas (jOmukai et al.ll2008l ). Our results on chemical enrichment suggests that quasi-stars 
are likely to be formed only in a redshift window between z ~ 20 and z ~ 10. At the lower 
end of this redshift window quasistars might only live in regions that lack primordial galaxies 
and thus have a relative suppression of structure formation compared to the average over 
the whole Universe. The environment where the progenitors of bright high redshift quasars 
live is unlikely to qualify as one of such regions. If we assume that a brig ht, rare, z = 6 



QSO is located at the center of a dark matter halo of mass M = 5 x 1O 12 M ( jSpringel et al. 



20051 ). then we derive that its first progenitor with T V i r > 10 4 K has been formed at z > 26 
with a confidence level greater than 0.9999. Such progenitor halo has a probability of being 
metal enriched greater than 98%, so it only has a small chance to host a quasi-star. Note 
however that quasistars can still lead to the formation of a bright z ~ 6 QSO if they are 
first formed in a relatively void region at z ~ 15 and then merge by z ~ 10 into the gas 
rich environment of the main QSO progenitor halo. This scenario allows the black hole seed 
from the quasistar to grow with maximal efficiency In fact, the BH starts with a large 
mass from the initial direct collapse and then it enters a gas rich region where it can grow 
continuously at near Eddington limit. As a z = 6 bright QSO halo consists of material 
originating in a sphere with comoving radius larger than 3Mpc, this is not unlikely. In fact 
a wind traveling at speeds of 30km/s only covers a fraction of this distance in half a billion 
years. However a detailed modeling can be obtained only through the study of QSO merger 
trees that include information on the spatial distribution of the progenitors, such as those 
given by cosmological simulations. 



5. Population III stars and the Reionization of the Universe 



The role of Population III stars in the reionization of the Universe has been much 
debated in the last years, especially after the first year WMAP data releas e, which included 
a high optical dept h to reionization r P = 0.17 ± 0.03 (ISpergel et al.l 120031) . Several models 
had been proposed (Venkatesan et al.ll2003l ; ICenll2003l ; IWyithe &: Loebll2003l ; iHui fc Haiman 



20031 ; IStiavelli et al.l 12004 ) and many of these included a significant contribution from first 
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stars to r e . With the latest WMAP data release, the optical depth to reionization is rather 



low (t p = 0.084 ± .016 implying an instantaneous reionization redshift z r 



10.0; see 



Komatsu et al. 2008) and its major contribution comes from complete ionization after z = 6 
(jShull &: Venkatesan!l2008l ). The contribution from higher redshift is lim ited to Ar R = 0.03 ± 
0.02, providing an upper limit to the luminosity of primordial galaxies (jShull fc Venkatesan 



2008|). 



From our study it appears that the contribution of Population III stars to the total 
budget of reionizing photons is limited if only one star per halo is formed, even neglecting 
negative radiative feedback. Despite the fact that first stars are more than one order of 
magni tude more efficient at producing ionizin g; photons per unit mass than Population II 
stars (ITumlinson &: Shulll l2000l ; ISchaererl 120031 ) their overall star formation rate in our stan- 
dard model is significantly lower for z < 20. Based on our standard model and assuming a 
Population III formation rate of lO _5 M Mpc" 3 yr~ 1 from z = 35 to z = 10, we obtain that 
about n^Ax 10 65 Mpc" 3 ionizing photons are emitted by metal free stars. This falls short of 
the number density of hydrogen atoms Uh ~ 7 x 10 66 Mpc~ 3 . Thus after taking into account 
the effect of dumpiness of the IGM and of recombination, it is clear that Population III stars 
can only reionize a minor fraction of the hydrogen atoms, even if the escape fraction is near 
unity. A large escape fraction is in fact possible in minihalos (M < 10 6 M Q ), but is likely 
significantly smaller in larger halos (M < 1O 6 M ), wh ere the H II region may remain con- 



fined well within the virial radius of the host halo (see IWhalen et al.l 12004 ; iKitayama et al. 



2004T ). The number of ionizing photons produced is smaller than the number of hydrogen 



atoms even in our model with no negative feedback (see Fig. [3}. In order for Population 
III stars to be a significant agent of reionization multiple PopIII stars must be formed in 
a single halo. For our model with ep op ni = 0.005, we obtain a cumulative ionizing photon 
production of ~ 4 x 10 67 down to z = 10. Such number of photons starts to become sufficient 
to contribute to reionization even for a relatively low escape fraction (f esc ~ 0.1). Certainly 
Population III stars are major agents of reionization if their star formation efficiency goes up 
one order of magnitude to ep op ni = 0.05 (see Fig. [7]). Note that in both these scenarios with 
multiple metal free stars per halo, the main sources of reionization are primordial galaxies in 
halos with T vir > 10 4 K. In fact, the reionizing photon budget from Population III remains 
significant even when only one star per minih alo is formed, but cluste rs of metal free stars 
are allowed in larger halos (see Fig. [HJ see also lHaiman fc Bryan! 120061 ). 
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6. Conclusions 

In this paper we present a model for the star formation rate of metal free (Population III) 
and second generation (Population II) stars during the Dark Ages of the Universe, at z > 10. 
The model relies on dark matter halo mass function coupled with analytical prescription for 
cooling and collapse of gas clouds. Our model includes radiative Lyman Werner feedback, 
which can suppress star formation in minihalos, and self-enrichment feedback, which marks 
the transition from metal free to Population II stars. 

Thanks to our novel treatment of chemical enrichment, based on the formalism devel- 



oped in lTrenti fc Stiavellil (120071 ). we show that halos with a virial temperature T V i r > 10 4 K 
are most likely to host a second generation of stars, formed from gas enriched to a metal- 
licity Z > 1O _4 Z by a progenitor Population III star in a minihalo at a higher redshift. 
Metal free stars can form in halos with T v i r > 10 4 K only once the cooling of gas in mini- 
halos is strongly suppressed by radiative Lyman Werner feedback, which in our reference 
model happens at z < 20. If only one Population III star forms per dark matter halo, then 
their number is dominated by those formed in minihalos with a peak star formation rate of 
pa 10~ 5 M Q Mpc _3 yr~ 1 at z pa 20. This prediction is robust and does not depend on the detail 
of the model. In fact, the negative radiative Lyman Werner feedback acts as a self-regulator 
of star formation in minihalos keeping variations of the star formation rate limited when the 
feedback efficiency or the halo mass function is changed. 

The metal enrichment also leads to an early rise of the star formation rate of Population 
II stars. By redshift z < 26 their SFR is higher than that of Population III stars and steadily 
rises as the redshift decreases. In our model we do not include positive radiative f eedback 
that can promote H 2 formation in the neighborhood of a first star (jRicotti et al. 2001 ). If this 



is the case, then the transition to metal enriched stars in halos with T V i r > 10 4 K is expected 
to be even more solid, because multiple Population III stars in clustered minihalos can 
pollute to a higher metallicity the gas that later constitutes a T vir > 10 4 K halo. The metal 
enrichment probability from minihalo pollution decreases once a strong LW background is 
in place, so at redshift z < 15 halos containing pristine gas with a mass M > 2 x 10 7 M Q are 
possible, provided that winds from protogalaxies, absent in our model, are not too efficient in 
polluting the I GM. Redshift z pa 15 m ight thus be the most favorable period for the formation 



of quasi-stars (jBegelman et al.ll2008l ). which are however expected to reside preferentially in 



underdense environments, where winds are more unlikely to be present and pollute the IGM. 

For this standard scenario the contribution to reionization given by Population III stars 
is only indirect (they enrich the IGM and allow Population II to form). In fact the cu- 
mulative number of ionizing photons they produce falls shorter of the number of hydrogen 
atoms. Metal-enriched stars in the first galaxies are thus expected to be the main agents 
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qf reionization, even though the rapid decrease of the galaxy luminosity function at z > 6 



( IBouwens fc Illingworthll2006l ; lOesch et al.ll2008l ) casts some doubts on this scenario 



A main change in the Population III star formation rate, which increases their contri- 
bution to reionization, can be introduced if one releases the assumption that only a single 
metal free star is formed per halo: by converting a fixed fraction of gas into primordial stars 
in minihalos the feedback mechanism is less efficient because as the critical gas mass needed 
for cooling is increases so is the number of stars produced per halo. Therefore in this scenario 
there is a constant growth of the Population III star formation rate in minihalos until this 
formation channel is suddenly inhibited because the minimum mass required to self-shield 
molecular hydrogen in the halo corresponds to a virial temperature above 10 4 K (see fig. |6]). 
Under this scenario Population III stars can easily produce a significant amount of ionizing 
photons and could account for a significant fraction of the optical depth to reionization orig- 
inating from z > 7. One issue to be addressed if multiple Population III stars are formed in 
a single halo is however the impact of local radiative feedback. The first massive star formed 
in a minihalo emits enough energy to completely photo-dissociate the gas in the halo, thus 
multiple stars can be formed only if there is a single star formation burst of limited time 
duration and very high efficiency This is a crucial issue that can be properly addressed only 
once cosmological simulations of Population III star formation will be able to go past the 
formation of the fi rst protostellar core and follow multiple episodes of star formation (see 



Wise fc Abelll2008l for promising progress in this direction). 
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Fig. 1. — Star Formation during the reionization epoch predicted by our standard model for 
Population III stars. Upper left panel: star formation rate versus redshift for Population 
III stars in minihalos (solid black line), Population III stars in more massive halos, with 
T V i r > 10 4 K (long dashed blue line) and for Population II stars formed out of metal enriched 
gas (short dashed red line). Upper right panel: minimum dark matter halo mass required to 
form a Population III star via H 2 cooling (solid black line) and in halos with T vir > 10 4 K 
(long dashed blue line) in function of redshift. Lower left panel: J21 flux in function of redshift 
with contribution from Population III stars (black shaded area) and from Population II stars 
(red shaded area). For z > 22 Population III stars are the main source of radiative feedback. 
Lower right panel: Probability of metal enrichment via progenitor polluti on for a halo with 



T vir = 10 4 K. The results have been obtained using a WMAP5 cosmology, a lSheth &: Tormen 



(119991 ) halo mass function and our model includes cooling and J21 feedback. The star 
formation rate is one star per halo for Population III stars and 5 x 10~ 3 for Population II 
stars. 
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Fig. 2. — Like in Fig. [T] but for our model with a reduced escape fraction, which implies 
a less efficient LW feedback. In this case the formation rate of H2 Population III stars is 
increased, while Population III stars in halos with T V i r > 10 4 i^ are partially quenched by 
chemical enrichment compared to fig. [TJ 
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Fig. 3. — Like in Fig. [T] but for our model with no LW feedback. Note that in this model 
J21 does not influence star formation, hence it is not shown. In this model quenching of 
Population III stars in the more massive T„j r > 10 A K halos is further enhanced compared to 
Fig.d 
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Fig. 4. — Like in Fig. [T] but for our model with a strong external radiative field (see 
eq. [HO This model yields similar results to our standard model in Fig. [T] because once the 
LW feedback is above the critical threshold necessary to quench star formation in minihalos 
its further increase has a modest additional effect. 
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Fig. 5. — Like in F i g, [ph ut for our model where the halo mass function is computed using the 
Press &: Schechterl ( 119741 ) formula. The star formation is suppressed at very high redsh i ft, but 
by z ~ 30 the model closely resembles our standard model where the lSheth fc Tormenl (119991 ) 
mass function is used. In fact, negative feedback acts as a self- regulator of the Population 
III star formation rate. 
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Fig. 6. — Like in Fig. [T] but for our model with multiple Population III stars allowed to form 
in the same halo (with efficiency ep opU j = 0.005). 
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Fig. 7. — Like in Fig. [T] but for our model with multiple Population III stars allowed to form 
in the same halo and have an enhanced efficiency (epopin = 0.05). 
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Fig. 8. — Like in Fig. [T] but for our model with multiple Population III stars allowed to form 
in halos with T vir > 10 4 K (with efficiency ep<ypiii_T vir =\§ A K = 0.005) and a single star per 
minihalo. 
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Table 1: Main Parameters for Star Formation and Feedback 



Model 


Halo MF 


tPopIII 


tPopII 


fesc 


Standard 


ST99 


1/halo 


5 x 1(T 3 


0.5 


LowEscape 


ST99 


1/halo 


5 x 1(T 3 


0.05 


NoEscape 


ST99 


1/halo 


5 x 10~ 3 


0.0 


Ext_J21 


ST99 


1/halo 


5 x 1(T 3 


0.5 


PS 


PS76 


1/halo 


5 x 10~ 3 


0.5 


MultiPopIII 


ST99 


5 x 1(T 3 


5 x 1(T 3 


0.5 


MultiPopIII_high_eff 


ST99 


5 x 1(T 2 


5 x 1(T 2 


0.5 


MultiPopIILLy-a 


ST99 


H 2 : 1/halo; T > lO 4 ^: 5 x 10~ 3 


5 x 10~ 3 


0.5 



Note. — Summary of the parameters that we modify from run to run in our star formation model. The 
first column identifies the model, the second reports the halo mass function used (ST99: Sheth & Tormen, 
PS76: Press & Schechter). The third and the fourth columns contains the star formation efficiency for 
Population 111 and for Population II stars. The last column the escape fraction. Model Ext_J21 includes an 
external radiative field (see Eq. [TO)) . 



